Shear flow pumping in open microfluidic systems 
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We propose to drive open microfluidic systems by shear in a covering fluid layer, e.g., oil covering water-filled 
chemical channels. The advantages as compared to other means of pumping are simpler forcing and prevention 
of evaporation of volatile components. We calculate the expected throughput for straight channels and show that 
devices can be built with off-the-shelf technology. Molecular dynamics simulations suggest that this concept is 
scalable down to the nanoscale. 



Progressive miniaturization and integration of chemical 
and physical processes into microfluidic systems, so-called 
"chemical chips", is expected to permit cheap mass produc- 
tion and to facilitate the handling of much smaller quantities 
than standard laboratory equipment 01 0] ■ Most currently 
built microfluidic systems contain the liquids in closed pipes 
and use body forces (e.g., gravity or centrifugal forces) |4|], 
electroosmosis, capillary wicking, or simply applied pressure 
to generate flow. Driving these systems is relatively straight- 
forward, but there are major drawbacks. Small pipes get eas- 
ily clogged (in particular while processing biological fluids 
which contain large molecules such as proteins or DNA), driv- 
ing gets increasingly difficult with reduced cross section, and 
fabrication more challenging. 

Open microfluidic systems are a second line of develop- 
ment which can overcome these problems. As demonstrated 
experimentally Q, 0] and theoretically 0], fluids can be 
guided on flat solid surfaces by wettability patterns, e.g., hy- 
drophilic stripes on an otherwise hydrophobic substrate. Sev- 
eral methods to generate flow in such "chemical channels" 
have been discussed: capillary forces (i.e., wicking into chan- 
nels or motion due to wettability gradients) |5], electrowet- 
ting 0, [l(ill . thermally generated Marangoni forces (e.g., in 
arrays of local heaters) II ll 1 1 211 - body forces, surface acoustic 
waves flUl . or combinations of these. While capillary forces 
can only induce flow over relatively short distances, electrode 
and heater arrays require complex sample structures. Surface 
acoustic waves can only drive drops which are large com- 
pared with the wavelength (in general on the micron scale), 
and high rotation velocities are required to obtain sufficiently 
strong centrifugal forces. 

As illustrated in Fig.^we propose to use shear, generated in 
a covering fluid layer (e.g., oil on water), to induce flow over 
long length scales, e.g., across the whole device. A top plate 
moving with respect to the bottom patterned substrate gener- 
ates shear in the oil which fills the gap between the two plates 
and, in turn, the oil will induce flow in a chemical channel 
covered with, e.g., water. This principle can be realized in a 
device with two circular plates rotating relative to each other 
and the chemical patterning can be implemented with stan- 
dard printing techniques Il4 Il51 . Design and scalability of 



FIG. 1 : (a) Flow in a water filled chemical channel (dark stripe on the 
static plate) can be induced by shear in the oil filling the gap between 
the static bottom plate and the moving top cover, (b) In a device the 
motion of the top with respect to the bottom plate can be realized 
by two circular discs rotating relative to each other. The pattern on 
the bottom plate indicates two reservoirs (circles) connected by two 
channels with two interconnections. 



shear driven systems is similar to body force driven devices 
and the same type of surface tension driven pearling instabil- 
ities occur in straight channels |8]. In both cases pearl for- 
mation increases throughput through the channels. This driv- 
ing principle can be extended easily to water confined in open 
topographic channels UM, although it will be more efficient 
in open chemical channel systems as there the water has less 
contact with the substrate (due to the lack of sidewalls). 

In the following we calculate the throughput for a shear 
driven straight chemical channel in two simple cases. Us- 
ing molecular dynamics simulations of a system of two im- 
miscible model fluids we demonstrate that shear driving can 
be downscaled to the nanometer scale. We conclude by dis- 
cussing how the principle of shear driving can be realized in 
a device of the size of a miniature computer hard disk (micro- 
drive). 

In order to estimate the possible throughput we consider the 
simple case of a straight and homogeneously filled chemical 
channel aligned parallel to the shear flow in the covering oil 
layer. We assume that the three phase contact line is pinned 
at the channel edges. In microfluidic devices the Reynolds 
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number is typically small so that one has Stokes flow. For 
two situations one can solve the corresponding hydrodynamic 
equations analytically: (i) the water-oil interface meets the 
substrate surface at a contact angle of 9 = 90°, (ii) the vis- 
cosities of the oil and water are equal. 

Figure |2 shows a cross section in a plane orthogonal to the 
chemical channel of width 2 a. Far from the channel we as- 
sume a constant shear rate 70. We seek solutions of the Stokes 
equation with the flow parallel to the channel (i.e., normal to 
the plane depicted in Fig.[2j. In this case, the pressure in the 
system is constant and the velocity component w parallel to 
the channel satisfies a Laplace equation 
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in both the water and the oil, with viscosity r\ w and r\ , re- 
spectively. By symmetry, the other velocity components are 
zero. For Couette flow in a homogeneous medium between 
two parallel moving plates (Fig. [Ha)) the solution of Eq. is 
w = 70 y, with a spatially and temporally constant shear rate 
70. In order to be able to use polar coordinates with the origin 
at the center of the chemical channel we take this solution as 
a boundary condition on an outer circle with radius R ;2> a, 
i.e., w\ r= ji = 70 R sin (Fig.|2j. We apply a no-slip bound- 
ary condition at the liquid-substrate interface for both fluids, 
i.e., ui|0 = o.tt = 0. At the interface between the two liquids 
the difference between the normal components of the stress is 
given by the Laplace pressure. For the geometry considered 
here this reduces to the condition that the interface has a con- 
stant curvature, i.e., it is a semicircle of radius r = a. The 
velocity field, w\ r=a -a = w\ r=a +o, as well as the tangential 
stress components, r\ w |^| r=a _ = r\ §^| r=a+0 , are contin- 
uous. For R/a 3> 1, the corresponding solution of Eq. Q is 



w(r, (f>) = 70 



< r < a 
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For R/a ^> 1 the flux (or throughput) Q = J A w d 2 r, with 
the channel cross section A, is given by 
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see Fig. |3 It increases with r\ jr\ w . The throughput can be 
characterized also by the apex velocity w*, i.e., the fluid ve- 
locity at <p = and r = a: 
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For r\ Q = rj w the flow field is not influenced by the presence 
of the liquid-liquid interface as long as the former is transla- 
tionally invariant along the channel. Accordingly, the velocity 
field is 



w{y) = 70 v- 




FIG. 2: Cross section of a chemical channel of width 2 a filled up to 
a contact angle 6 (here 6 = 90°) and filling height ho. In order to be 
able to use polar coordinates P = (x = r cos (j>, y = r sin </>) we 
fix the shear rate at r = R and consider the limit R/a — > 00. 
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FIG. 3: Normalized flux Q for 9 — 90° as a function of rj /r] w 
(full line) and for ?y D = r\ w as a function of ho / a (dashed and dotted 
line). The asymptotics for small and large arguments are indicated. 
For ho j a > 1 the liquid ridge is unstable with respect to pearling 
(dashed line). 



Eq. 10 implies as water flux through a channel with a cross 
section forming a circular segment with filling height ho(9) = 
a (1 — cos 0) j sin 9 (see Fig. [5} 
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so that Q vv (h = a) = Q 90 o (i] = r) w ) (see Fig.0. 

In order to suggest that shear driven microfluidics can be 
miniaturized to the nano-scale, we have performed molecu- 
lar dynamics simulations as a well suited approach in this 
regime. To a large extent we have used standard tech- 
niques with the implementation of the dynam- 
ics on wetting stripes following Ref. Ja]. The model sys- 
tem consists of two immiscible monatomic fluids and two 
types of wall atoms, with Lennard-Jones interactions Vij (r) = 
(5) 4e [(r/cr)~ 12 - (r/a)~ 6 ] truncated at r = 2.5 a, where r 
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is the interatomic separation, a measures the size of the repul- 
sive core, e is the strength of the potential, and Ay = A,, is a 
dimensionless parameter that controls the attraction between 
atoms of the four atomic species i and j. The temperature 
is set to T = l.Oe/ks by a Nose-Hoover thermostat. The 
liquids are confined between solid walls at a distance 50.7 c, 
each composed of a layer of fee unit cells (with lattice con- 
stant 1.55 a) whose atoms are tethered to lattice sites with a 
stiff linear spring. Periodic boundary conditions are applied in 
the two lateral directions. Couette flow is achieved by trans- 
lating the tether sites of the upper wall at a constant veloc- 
ity. Most of the interaction strengths Aj have the standard 
value 1.0, but Ay = if i and j refer to different liquids, or 
the inner liquid (the "water") and the non-wetting part of the 
wall, or the outer liquid (the "oil") and the wetting part of the 
wall, or to wall atoms. With these choices, the two liquids 
are completely immiscible. Because the interactions among 
atoms of each of the two species are identical, they have equal 
density and viscosity. A number of different system sizes have 
been simulated, all giving similar results; the largest contained 
about 1.25 million atoms in total, with wetting stripe dimen- 
sions 16.7a x 547(7. 

After a short equilibration period the inner liquid occupies 
a roughly cylindrical region pinned on one side to the wet- 
ting stripe (see Fig. Ela)). We note that the liquid-liquid in- 
terface is much smoother than the previously studied liquid- 
vapor interfaces above wetting stripes 0]. If the upper solid 
wall is moved parallel to the stripe, a two-fluid Couette flow 
results. Velocity profiles in two different vertical slabs paral- 
lel to the yz-plane are shown in Fig.|4jb). The velocities are 
averaged over the slab width of 4 a and the full system length, 
and a time interval of 250 r, with r = a \/m/e and the atom 
mass m. On the sides of the simulation box, where only the 
outer liquid and the non-wetting solid are present, the profile 
(o) varies linearly aside from a weak bending near the walls. 
In the center of the simulation box the velocity (+) shows 
some additional deviation from linearity in the inner liquid 
region. There is liquid-liquid interfacial slip (as expected for 
the strongly immiscible liquids used here 12011 ') and some pin- 
ning behavior at the wetting stripe. In both cases, the velocity 
near both the stationary and moving wall exhibits a no-slip 
behavior, although extrapolations of the velocity profile from 
the interior of the channel would give rise to an apparent slip. 

Under these circumstances the anchoring of the inner liquid 
at the wetting stripe is very stable, and we have observed that 
the inner liquid remains atop the stripe even if the forcing oc- 
curs at oblique angles. As an extreme example, if the top plate 
moves in the x-direction with the same velocity as above, the 
snapshot of the atomic positions in Fig.|^c) shows the inner 
fluid to be only modestly displaced to the side. However, in 
general, the interfacial shape is not stable, because the con- 
figuration is subject to a surface tension driven instability as 
long as the channel is long enough and filled to more than a 
90° contact angle 01 . In that case, a periodic surface undu- 
lation develops, as shown in Fig.|4ld), which eventually leads 
to propagating liquid pearls. 
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FIG. 4: (a) Cross sectional view of the molecular configuration at 
early times before flow commences. The inner stripe of substrate 
atoms (full boxes) attracts the inner liquid atoms (full circles), while 
the outer substrate atoms (open boxes) attract the outer liquid atoms 
(open circles). Atoms located within a slab of thickness 2.5 a are 
shown, (b) Velocity profiles of the shear flow at the side of the sim- 
ulation box (o) and at the center of the channel (+), in slabs cutting 
through the outer liquid and through the liquid ridge, respectively. 
The upper solid wall is moved in ^-direction at a velocity 0.3 <t/t. 
(c) Shear at an oblique angle with a component orthogonal to the 
wetting stripe results in a tilted shape, but the inner liquid ("water") 
is still transported along the channel, (d) Chemical channels filled to 
more than 9 = 90° are also unstable with respect to pearl formation 
under shear (snapshot of the liquid-liquid interface at t = 2500 r). 



In order to discuss the applicability of our concept (see 
Fig. Q, we estimate the throughput in a potential device. 
One prominent application is handling biological fluids which 
contain DNA fragments with a radius of gyration of about 
100 nm. We therefore consider a channel width 2 a = 1 /im 
and a plate spacing d = 0.1 mm. These are dimensions 
already considered for chemical channels in slit pores l2lll . 
Modern miniature hard disks (so-called Microdrives® with 
a storage capacity of up to 8 GB) have a disk radius p 
(see Fig. [flb)) of about one centimeter and rotate at v = 
3600 rpm = 60 Hz (see, e.g., Hitachi Global Storage Tech- 
nologies http://www.hgst.com). The whole device has a size 
of 42.8 mm x 36.4 mm x 5.0 mm so that building on this tech- 
nology might lead to microfluidic devices of similar dimen- 
sions. For p = 1 cm and u> — 2 it v = 400 Hz the relative 
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velocity of the discs at their perimeter is v = pto = 4m/s 
and the shear rate is 70 = v/d = 4 x 10 s . Assuming 
90° contact angles at the channel edges and that the viscosi- 
ties of the two fluids are equal Eq. @ renders an apex velocity 
u>9 o = 070 = 2cm/s. 

In an inertia-driven situation the velocity at the surface of 
a planar film of thickness a driven parallel to the substrate 
surface can be used as an order of magnitude estimate for 
the apex velocity w* = a 2 g w g/(2r] w ) with the mass den- 
sity g w of the inner fluid and the acceleration g. Thus for 
r\ w = 10~ 3 kg/(ms) and g w = 10 3 kg/m 3 an apex velocity 
w* = 2 cm/s requires an acceleration of g = 1.6 x 10 5 m s~ 2 . 
Such accelerations pto 2 can be achieved with centrifugal 
forces at a radius of 1 cm at 3.8 x 10 4 rpm as compared with 
a spin coater with 10 4 rpm, a turbo pump with 6 x 10 4 rpm, 
and a dentist's drill with 3 x 10 5 rpm. 

In order to assure laminar flow, the Reynolds number has to 
be small, both in the channel and in the oil. With the values 
for water one has Re u , = g w w* a/r) w = 1CP 2 in the channel. 
The relevant length scale for the oil is d and the velocity scale 
is p u) so that for g a as g w and r/ as lO 3 ^, one obtains a 
similarly low Reynolds number for the outer liquid, leading to 
laminar flow. 

Viscous heating of the oil is a possible concern. For the 
parameters discussed above, the power input into the device 
is of the order of r) w 7q dir p 2 /A as 13 mW so that cooling of 
the device is not expected to be a problem. 

In summary, we propose to drive flow in open microfluidic 
systems by shear in a covering liquid of equal or higher vis- 
cosity, e.g., water filled channels covered with oil. As a pos- 
itive side effect, the oil prevents evaporation. The analysis of 
the Stokes flow for channels filled to a 90° contact angle at the 
channel edges shows (see Fig.|3j how increasing the viscosity 
of the covering liquid leads to an enhanced channel through- 
put for the same shear rate in the covering liquid. Fig.[3]also 
shows how, in the case that the covering liquid and the liquid 
in the channel have the same viscosity, the efficiency of the 
shear driving increases with increasing filling level. This ob- 
servation is in agreement with the corresponding one in body 
force driven systems 0] . Our MD simulations have encour- 
agingly explored the potential to miniaturize the concept of 
shear driving down to the nano-scale. 

Typical shear driven microfluidic devices might consist of 
two small discs rotating relative to each other. One or both of 
the discs can have a pattern of hydrophilic stripes on their hy- 
drophobic surface. We discus how with available technology 
such a device can be built with the size of a miniature hard 
disc, i.e., a few centimeters wide and long and a few millime- 
ters thick. 

As expected, for a contact angle larger than 90° we ob- 
serve a surface tension driven pearling instability similar to 
the Rayleigh-Plateau instability as discussed in detail for body 
force driven systems (8fl. In MD simulations of a liquid com- 
pletely wetting (6> eq = 0°) a strip embedded into a non-wetting 
substrate ((9 eq = 180°) the liquid on the channel did not de- 
tach from the channel even if the shear was not aligned with 



the channel direction. The issue of the shear rate at which 
the liquid detaches from the channel is related to studies of 
moving droplets with contact angle hysteresis I22II23I1 . 

Our analysis suggests a significant application potential for 
shear driven open microfluidic systems, accompanied by fur- 
ther theoretical and experimental research. 
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